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Abstract 

We investigate the dependence of the Yang-Mills wave functional in Coulomb gauge 
on the Faddeev-Popov determinant. We use a Gaussian wave functional multiplied by 
an arbitrary power of the Faddeev-Popov determinant. We show, that within the resum- 
mation of one-loop diagrams the stationary vacuum energy is independent of the power 
of the Faddeev-Popov determinant and, furthermore, the wave functional becomes field- 
independent in the infrared, describing a stochastic vacuum. Our investigations show, 
that the infrared limit is rather robust against details of the variational ansatze for the 
Yang-Mills wave functional. The infrared limit is exclusively determined by the divergence 
of the Faddeev-Popov determinant at the Gribov horizon. 
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1 Introduction 



Recently Yang-Mills theory in Coulomb gauge has become the subject of intensive studies both 
on the lattice, refs. [I], [2] and in the continuum, ref. [3], [I], [7j, jHj. The Coulomb gauge 
is a physical gauge and in this gauge confinement is realized by the statistical dominance of 
the field configurations near the Gribov horizon, which gives rise to an infrared enhanced static 
color charge potential. 

For the calculation of static properties of continuum Yang-Mills theory, the Schrodinger equa- 
tion approach j^j seems to be most convenient. In refs. jHj, [Z], [E| the Yang-Mills Schrodinger 
equation was approximately solved in Coulomb gauge, using the variational principle. Differ- 
ent ansatze for the vacuum wave functional and different renormalization conditions have been 
used, and different infrared behaviours of the gluon and ghost propagators were obtained. One 
might wonder, whether the different results are a consequence of the different ansatzes for the 
wave functional. To answer this question, in this paper we consider a more general class of 
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wave functionals, which includes, in particular, the wave functionals previously used in refs. 
0, [7] and jH]. We will show, that the different ansatze used so far, have to yield the same 
unique infrared behaviour of the vacuum wave functional (at least to the order considered). We 
will also show, that in the infrared the wave functional becomes field independent, describing 
a stochastic vacuum. Furthermore, the infrared limit of the wave functional agrees with the 
exact vacuum wave functional in D = 1 + 1. 



2 The variational ansatz 



For the Yang-Mills vacuum we consider trial wave functionals of the form 

VIA 1 ] = j-^m^] , (1) 

where 

is the Faddeev-Popov determinant, which for later convenience has been normalized to J[A ± = 
0] = 1. Here D = d + A 1 - is the covariant derivative and A 1 - = A ±a T a denotes the gauge field 
in the adjoint representation. Furthermore (f)[A ] is a Gaussian wave functional defined by 

ftA 1 ] = Afexpi-SIA 1 ]) 

SiA^ = l -j 'd?x j 'dV4 ±a (x)<(x,xXV) (3) 

with M = jV(q!, uj) being a normalization constant to ensure (ip \ ip) = 1. The Faddeev-Popov 
determinant arises as Jacobian in the transformation from "cartesian" coordinates A°(x) to 
the "curvilinear coordinates" A^- a (x) satisfying the Coulomb gauge diAf- = 0, and defines the 
metric in the space of transversal gauge orbits A^-(x). Accordingly the scalar product in the 
space of transversal gauge orbits is defined by 

DA ± J[A ± ]V*lA ± }<S>lA ± } , (4) 

where the integration should in principal be restricted to the fundamental modular region jl], 
[Tij . The choice o;^ b (x, x') = 5ijS ab uj(x, x') and a = was used in refs. [7j, while a — | was 
chosen in ref. jSj. From eq. ([TJ it is seen, that for the latter choice (^[^4^] represents just the 
"radial" wave functional 1 . 

We wish to study the dependence of the vacuum Yang-Mills wave functional (1) on the power 
of the Faddeev-Popov determinant a, which can take, in principle, any real value as long as 

1 For a point particle in a s-state the wave function is of the form \l/(r) = where 4>{r) is the radial wave 
function and the Jacobian is given by J = r 2 
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^[A- 1 ] is normalizable. The integral kernel u as well as the parameter a have to be determined 
by minimizing the expectation value of the energy 

(H) = J DA ± J[A ± ]^*[A ± ]H^[A ± ] . (5) 



3 Minimization of the energy 

Inserting the explicit form of the wave functional eq. (0) into eq. (0) variation of the energy 
with respect to the kernel uj yields the "gap equation" : 

^-'(EK({i.-})-o. 

Here the first term arises from the variation of the normalization constant (5J\f/ 5uj = J\f(5S/ Su)) 
and {, } denotes the anti-commutator. Minimization of the energy (H) (jHJ) with respect to the 
power a yields the condition 



d(H) 
da 



2<lnJ[A J -])<lf) - ({in M^H}) = , (7) 



where we have used dN /da = Af(\n JfA- 1 ]). 

Consider now the structure of the Faddeev-Popov determinant (2), which obviously satisfies 
In JlA 1 - = 0] = 0. Furthermore, by definition (J2J) we have 

^ = -i>(<ar<»)>. w 

where 

G = {-DAY 1 (9) 

is the inverse Faddeev-Popov operator and T° k ' a (x.) = 5G~ X /5A^ a (x.) is the bare ghost gluon 
vertex [S]. Since r 0,a ~ T a (group generator in the adjoint representation) and T a occurs in G 
only in the combination A 1 - = A ±a T a it is clear, that the quantity (jHJ) has to be proportional 
to A- 1 , since trT a = 0. Therefore we find the representation 

In J[A X ] = J d 3 a;d 3 x'Cf/[A ± ](x,x / )A ±a (x)^f(x , ) (10) 

with some, not explicitly known functional C?- [A ]. In one-loop approximation the expectation 
value of eq. (fTUj) is given by 



(In JlA 1 ]) ~ J rf 3 xd 3 x'(Cf/[A ± ](x,x') 

[Ai a (x)Aj b (x')) . (11) 
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Furthermore, to this order we can neglect terms of the form (^j^A^ and ( sa^a 1 - ) anc ^ ^ nc ^ 
from for the curvature in orbit space [Hj 

> (xx1 _ 1/ ^ J \ 
' j 2\M>(x)<L4| 6 (x')/ 

= -(C^](x,x')) . (12) 
To this order we can also replace C[A] in (|TT|) by its expectation value (— x) yielding 

In J[A L ] = - J d 3 xd 3 x' X %(x, x') A^ a (x)Af b (x') . (13) 
Inserting eq. (|T5|) into eq. (JJJ) we find 

^ = "/ ^V X J(x,x')[2 (A ±a (x)Af (x')) <#> - ({A ±a (x)Af (x'),^})] (14) 
On the other hand for the Gaussian wave functional JI}, © we have 

V a (xKV), (15) 



^'f(x.x') > ' ' 



so that the equation (JHJ) becomes 

2 ^(x|xO = 2 (^(-KV)) <#> " ({A ±a (x)Af (x')F}) . (16) 
Comparison of eqs. (pPIj) and (fTB^) yields 



da 



-2/^V X -(x, X ' )s ^. (17) 

Thus stationarity of the energy with respect to ^ b (x) , 5(H) /Sou = implies also stationarity 
with respect to a,d(H)/da = 0. Let us emphasize, that eq. (JT7J) is exact to one-loop order in 
the equation of motion (i.e. to two- loop order in (H)). 



4 The energy functional 

The above obtained result can be also immediately infered from the explicit expression of the 
expectation value of the Yang-Mills Hamiltonian in the state which is given by 

(H) = E k + E B + E c (18) 
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where d(k) and /(k) are the ghost and Coulomb form factors defined in ref. [H] and x is the 
scalar curvature defined in terms of the curvature tensor (J 12)) by 



ub i 



(19) 



X(k) = ^ / 7^ f 1 - (kq) 2 l ' (k ' q) 1 q) (20) 
Ay ' 4 J (2tt) 3 L v h; J (k- q ) 2 v ; 

with tfci(x) = <5fcz — d^di/d 2 being the transversal projector 2 . Furthermore 

fl(k) = w(k) - (2a - l) X (k) (21) 

is the inverse of the gluon propagator 

(A ±a (k)Af (-k)> = V^O-^k) . (22) 

Note, that the curvature x(k) (J2DJ1 is entirely determined by the ghost form factor d(k) and does 
not depend on u;(k). The energy ()18|) depends on a and ct>(k) only through the combination 
Q(k) = uj(k) — (2a — l)x(k). From this fact immediately follows, that eq. © implies eq. (0), 
so we find again, that the wave functional (0) which minimizes the energy is independent of a. 
In fact, since (H) (|18|) depends on u and a only through the combination Q = uo — (2a — l)x 
it suffice to minimize the energy with respect to Q. The resulting gap equation 3 also depends 
only on Q and its solution is independent of a. This shows, that the infrared behaviour of the 
gluon propagator (A- 1 A- 1 ) ()2*2*|) is independent of the power a of the Faddeev-Popov determinant 
assumed in the wave functional (0. Therefore we are free to choose a for our convenience, for 
example, a — |. This choice has the technical advantage, that f2(k) = cu(k), which allows a 
straightforward application of Wick's theorem in the calculation of expectation values. 
In this context let us also mention, that the choice a = | in eq. (0) yields the wave function 
used by the present authors in ref. [S] , while the wave function used in refs. P , [I] corresponds 
to the choice a = 0. Inspite of the different wave functions chosen in refs. jS], [Zj and ref. jHj 
the same infrared behaviour of the gluon propagator should be obtained in one-loop order, as 

2 The ghost and Coulomb form factor, d(k) and f(k), satisfy the Schwinger-Dyson equation derived in ref. 
jS] with u> replaced by fl 

3 Although the curvature (121 )[) does not explicitly depend on uj it depends implicitly on u> via the ghost 
form factor d(k). However, also the ghost form factor d(k) depends on u> only through il. This dependence is, 
however, a higher order effect and to one-loop order |^ or & can be neglected. Then the gap equation = 
is exactly the one obtained in ref. [S] with u> replaced by f2 . 
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shown above, provided the same renormalization condition is used. However, while refs. |SJ, [7j 
finds an infrared finite gluon propagator, we find an infrared vanishing gluon propagator jS]. 
Two sources of the different behaviours obtained in ref . jH] , and refs. jB] , [7] come to mind: i) 
different choices of the renormalization condition and ii) different treatments of the curvature 
of orbit space. In ref. jH] we choose the socalled horizon condition ^3] 

cT\k -> 0) = (23) 

as renormalization condition while refs. jE], [Zj require the kernel u in the Gaussian wave 
functional to be infrared finite u(k — > 0) = const. 4 . Furthermore, while the curvature of orbit 
space x was fully included in ref. jHJ , it was completely neglected in j^] and ignored 5 in the 
Coulomb energy in the numerical calculations of ref. [7_j. As will be shown in the next section 
ignoring the curvature in the Coulomb energy will change the infrared behaviour of the wave 
functional. It was already observed in ref. [HJ, that the full inclusion of the curvature is vital 
for the infrared limit of the theory. This is consistent with the observation in Landau gauge, 
that in the Schwinger-Dyson equations the ghost loop is by far more important than the gluon 
loop !■ 



5 The Yang-Mills wave functional in the infrared 

With the relation (fT"3*|) the wave functional (0) becomes 

^[A^ ~ e <* f ^xA^U A± " A± . (24) 

Furthermore, the solution of the gap equation eq. (jUJ) is such, that in the infrared 

x (k -> 0) = Q(k -> 0) , 2ax(k -> 0) = cu(k -> 0) (25) 

holds. This is an extension of the relation x(k — >• 0) = uj(k — > 0) found in ref. jH] for a = |. 
With the relation (|25jl the vacuum Yang-Mills wave functional becomes in the infrared 

= 1 . (26) 

In ref. pLOj this wave functional was assumed in the infrared regime, for sake of simplicity. We 
have thus shown, that, to one-loop order, eq. ()26|) is the correct wave function in the infrared. 
The infrared wave functional ^/[A- 1 ] = 1 means, that gauge fields at distant positions x, x', |x — 
x'| — > oo are completely uncorrelated. Accordingly, the gluon propagator (y4 J -(x)A- L (x')) has 
to vanish rapidly in the infrared |x — x'| — > oo, which is in agreement with the findings of 
ref. jBJ. Thus, the wave functional ^/[A- 1 ] = 1 describes a stochastic vacuum, in which color 
(correlation) cannot propagate over large distances. This is nothing but color confinement. In 

4 In D = 2 + 1 wc find a self-consistent solution to the coupled Schwinger-Dyson equations only when we 
impose the horizon condition J22J). 

5 In the formal part of |7] the curvature was fully included. 
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this sense, the infrared wave functional supports the picture of a stochastic Yang-Mills 
vacuum [TTj . 

Let us also stress, that in view of the relation (J25j) the infrared behaviour of the gluon propagator 
()22|) is exclusively determined by the curvature x in orbit space. Furthermore from eqs. 
()26|) and (jl3j) follows that in the infrared limit the vacuum expectation values are described by 
an Gaussian ensemble 

{■■■) = J VA ± J[A ± }--- = J DA ± ---e~I A±xA± . (27) 

Finally, let us also note that the relation ()25|) holds independent of the employed renormaliza- 
tion condition as long as the curvature x(^) is infrared divergent. Since the Faddeev-Popov 
determinant vanishes on the Gribov horizon, which contains the infrared dominant field con- 
figurations, from eq. (fTH|l follows that the curvature has, indeed, to be infrared divergent. The 
condition (|23j) is, however, lost when the curvature is neglected in the Coulomb energy as done 
in ref. jH], |Zj- Given the infrared singular behaviour of x{k 0) t ne condition ()25|) implies 
that for a 7^ (in particular for a = | [Hj) the variational kernal u{k) in the Gaussian ansatz 
(J2J) has to be infrared singular, while the choice a = 0, [Zj can tolerate an infrared finite 
u(k). 



6 Yang- Mills theory in D = 1 + 1 

Let us test the above result in 1 + 1 dimension, where Yang-Mills theory can be solved exactly 
on a torus and reduces to quantum mechanics in curved space. 

Implementing the Coulomb gauge d±Ai = 0, there is only a constant gauge field Ai(xi) = const, 
left, which can be diagonalized in color space by exploiting the residual global gauge freedom 
U, not fixed by diAi = 0. Defining the remaining quantum mechanical degree of freedom, a, 
by 

gA 1 L = gLAf-^ = U~tW , (28) 



where L is spatial extension of the torus the Faddeev-Popov determinant becomes 

J(a) = sin 2 a. (29) 

The Gribov horizon occurs at a = utc and the fundamental modular region is obviously given 
by < a < it. Furthermore the Yang-Mills Hamiltonian in the variable a is given by 

g 2 L 1 d . 2 d , , 

H kin = 3 — r^— -j- sin a— . (30) 

8 sin a da da 

In one spatial dimension there is no magnetic field and no dynamical gluon charge (— A1 b U1 = 
0), since the gauge field has only one (non-zero) color degree of freedom. Accordingly, the 
Coulomb term of the Yang-Mills Hamiltonian [3] vanishes in the absence of external color 
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charges. 

With the ansatz 



V k (a) = -^=(f) k {a) = -?—(t> k {a) , (31) 

sin a 



which corresponds to the choice a = | in eq. (JJJ), the Schrodinger equation H^ k = E k ^/ k 
reduces to 

- = (E k + ^] Ma) , (32) 



whose solution is given by 6 



a 2 L 

i fc (a) = sin(fca) , E k = y —{k 2 - 1) . (33) 



In the continuum limit L — > oo only the vacuum state k = 1 survives (E± = 0), while all excited 
states k > 1 aquire an infinite energy and are thus frozen out. The vacuum wave function is 
given by (jSU), (El 

*k=i(a) = l, (34) 

which is precisley the infrared limit of the vacuum Yang-Mills wave functional in D = 3 + 1 
found above (see eq. (|26))). Note also that the radial wave function <p k (a) ()33|) vanishes on 
the Gribov horizon a = ni\ to compensate for the vanishing of the Faddeev-Popov determinant 
J [a] (|2*nj) . just like in the D = 3 + 1 dimensional case where (for a — |), the radial wave 
functional (^[^l" 1 ] © vanishes in the infrared due to the infrared divergence of u(k — > 0). 



7 Summary and Conclusions 

We have studied the variational solution of the Yang-Mills Schodinger equation in Coulomb 
gauge for a class of wave functionals ((TJ) consisting of a Gaussian and an arbitrary power (—a) 
of the Faddeev-Popov determinant. We have found, that up to one-loop in the gap equation 
(i.e. two loops in the energy) the stationary solution is independent of this power a. The same 
is true for the transversal gluon propagator (|2*2*jl which is exclusively determined by the self- 
consistent solution n of the gap equation = 0. This solution is independent of the choice 
of a. Different choices of a will lead to different kernels u (with possibly different infrared 
behaviours) in the wave functional (|3*|). But this will not affect the gluon propagator (|22jl. 
Furthermore in the infrared the Yang-Mills vacuum wave-functional becomes field-independent 
describing a stochastic vacuum, in which color cannot propagate over large distances. The 
infrared limit of the wave functional becomes exact in D — 1 + 1. 

Our investigations show, that the infrared behaviour of Yang-Mills theory in Coulomb gauge 
is rather robust with respect to changes in the variational ansatze for the wave functional as 
long as the curvature in orbit space induced by the Faddeev-Popov determinant is properly 
included. 

6 This solution was previously found in ref. |13j 
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